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For  seal-less  type  solid  oxide  fuel  cells,  its  power  generation  characteristics  and  distribution  of  the  gas 
composition  depend  on  not  only  the  electrochemical  reaction,  but  also  complex  kinetics  and  transport 
phenomena,  because  the  internal  reforming  reaction  and  the  water-gas  shift  reaction  take  place  together 
with  reverse  diffusion  of  the  ambient  gas  from  the  surroundings  of  the  cell.  The  purpose  of  this  paper  is 
to  theoretically  explain  the  experimental  results  of  the  anodic  concentration  profile  of  gaseous  species 
previously  reported  in  a  practical  seal-less  disk-type  cell  which  used  pre-reforming  methane  with  steam 
as  a  fuel.  A  numerical  model  that  takes  into  account  the  transport  phenomena  of  the  gaseous  species 
and  the  internal  reforming  reaction  with  the  water-gas  shift  reaction  together  with  the  assumption  of  the 
cell  outlet  boundary  condition  was  constructed  to  numerically  analyse  the  gas  composition  distribution 
and  power  generation  characteristics.  Numerical  analyses  by  the  model  were  conducted  for  the  several 
cases  reported  as  the  experiment.  The  calculated  results  in  the  anode  gas  concentration  profile  and  in  the 
voltage-current  characteristics  show  good  agreement  with  the  experimental  data  in  every  case,  and  then 
the  validity  of  the  simulation  model  was  verified.  Therefore,  the  model  is  useful  for  a  seal-less  disk-type 
cell  which  is  operated  by  a  fuel  including  non-reformed  methane. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  have  attracted  much  attention  due 
to  their  high  efficiency  and  cleanliness,  and  they  are  expected  to  be 
one  of  the  most  favorable  electrical  power  generation  technolo¬ 
gies  in  the  near  future.  Recently,  SOFC  systems  driven  by  town 
gas  have  been  intensively  developed  by  several  companies  using 
seal-less  cell  structure,  where  the  outlet  edges  of  the  gas  channels 
of  the  anode  and  cathode  are  open  to  the  ambient,  owing  to  its 
advantageous  features  in  gas  sealing  and  thermal  stress  problems 
[1-3]. 

In  the  SOFC  systems,  an  internal  reforming,  i.e.,  reforming  of  the 
fuel  in  the  cell,  is  often  employed  in  conjunction  with  pre-reforming 
due  to  its  high  operating  temperature.  Generally,  the  power  gen¬ 
eration  characteristics  of  cells  depend  on  the  gas  composition  on 
the  anode  side  and  cathode  side  of  the  cells.  In  the  seal-less  type 
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cells,  the  internal  reforming  reaction  and  the  water-gas  shift  reac¬ 
tion  take  place  together  with  reverse  diffusion  of  the  ambient  gas 
from  the  surroundings  of  the  cell,  i.e.,  diffusion  opposite  to  gas  flow 
direction.  Therefore,  the  distribution  of  the  gas  composition  in  the 
cell  and  power  generation  characteristics  depends  on  not  only  the 
electrochemical  reaction,  but  also  complex  kinetics.  Many  simula¬ 
tion  studies  have  been  conducted  to  aid  in  development  of  SOFCs 
[4-19]. 

To  the  best  of  our  knowledge,  however,  very  few  numerical 
results  have  been  reported  on  the  power  generation  characteristics 
compared  to  the  measurement  of  the  gas  composition  distributions 
because  the  concentration  profile  of  the  gaseous  species  in  a  cell  has 
been  seldom  measured.  For  the  seal-less  disk-type  cell  using  pure 
hydrogen  as  a  fuel,  it  has  been  reported  by  Momma  et  al.  [17]  that 
the  influence  of  the  reverse  diffusion  of  the  ambient  gas  is  clari¬ 
fied  based  on  both  the  experimental  and  theoretical  analyses,  the 
results  of  which  have  shown  good  agreement.  On  the  other  hand, 
there  has  never  been  such  practicably  useful  report  on  the  cell  using 
reformed  gas  as  the  fuel. 

In  one  of  the  recent  papers  dealing  with  seal-less  disk-type 
SOFCs,  Momma  et  al.  [20],  reported  the  experimental  results  for  the 
concentration  profile  of  gaseous  species  along  the  flow  direction  in 
the  anode  gas  channel  of  the  single  cells  which  were  performed 
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Nomenclature 

c  molar  density  (mol  m-3 ) 

cCo  CO  concentration  in  the  mixture  gas  (mol  m-3 ) 

Ch2o  H20  concentration  in  the  mixture  gas  (mol  m-3) 

Cj  mole  reaction  quantity  due  to  cell  reaction  for  each 

unit  time  and  unit  volume  (mol  s-1  m-3) 
cref  mole  reaction  quantity  due  to  reforming  reaction  for 
each  unit  time  and  unit  volume  (mol  s-1  m-3) 
cSf  mole  reaction  quantity  due  to  water-gas  shift 
reaction  for  each  unit  time  and  unit  volume 
(mols-1  m-3) 

Dab  mutual  diffusion  coefficient  of  species  A  and  B 

(m2  s-1) 

dan  thickness  of  anode  current  collector 

dca  thickness  of  cathode  current  collector 

De  effective  diffusion  coefficient  (m2  s-1 ) 

Dj  m  diffusion  coefficient  of  species  i  in  multi-component 
mixture  gas  (m2  s-1 ) 

Dy  diffusion  coefficient  of  the  pair  of  1  -j  species  in  a 
binary  mixture  (m2  s-1 ) 

Dim  diffusion  coefficient  when  element  1  diffuses  in  the 
mixture  gas  m 

Ear  activation  energy  (J  mol-1 ) 

A Earef  activation  energy  ( =57,840 J  mol-1 ) 

A EaSf  activation  energy  (=95,000J  mol-1 ) 

AG^2q  Gibbs  free  energy  change  in  reaction  H2  +  02  =  H20 
at  the  standard  pressure  (0.1  MPa)  and  reference 
state  (J  mol-1 ) 

F  Faraday  constant  (C  mol-1 ) 

Je  electronic  current  density  which  flows  in  electrolyte 

(Am-2) 

Ji  ionic  current  density  which  flows  in  electrolyte 

(Am-2) 

Jo  current  density  which  flows  to  an  external  circuit 

(Am-2) 

k  Boltzmann  constant  (=1.3807  x  10-23  J  K-1 ) 

I<Sf  equilibrium  constant  of  water-gas  shift  reaction 
Mi  molecular  weight  of  species  i  (kg  mol-1 ) 

P  pressure  (Pa) 

Pch4  partial  pressure  of  methane  in  the  mixture  gas 
Pi  partial  pressure  of  chemical  species  i 

R  gas  constant  (=8.314J  mol-1  K-1 ) 

rceu  radius  of  the  cell 

Rceu  cell  resistance  ( £2  m2 ) 

Rceii0  cell  resistance  at  1073  K  ( £2  m2 ) 

Re  ohmic  resistance  of  electrolyte  by  electronic  con¬ 

duction  (£2  m2) 

Ri  ohmic  resistance  of  electrolyte  by  ionic  conduction 

(£2  m2) 

Rr  reactive  resistance  of  the  electrodes  (£2  m2) 

Rs  series  resistance  of  the  electrodes  (£2  m2) 

rcH4  rate  of  the  reforming  reaction  (mol  s-1  kg^j.) 

rc o  rate  of  the  water-gas  shift  reaction  (mol  s_1  kg”^) 

T  temperature  (I<) 

Tceu  temperature  of  the  cell 

T'ref  outlet  gas  temperature  of  pre-reformer  (I<) 

Vcen  cell  operating  voltage  (V) 

u  mass  average  velocity  (m  s_1 ) 

Uf  fuel  utilization 


Greek  letters 

a ref  utilization  factor  by  which  the  catalyst  in  the  anode 

current  collector  actually  contributes  to  the  reform¬ 
ing  reaction 

aSf  utilization  factor  by  which  the  catalyst  in  the  anode 
current  collector  actually  contributes  to  the  water- 
gas  shift  reaction 

san  porosity  of  the  anode  current  collector 
sca  porosity  of  the  cathode  current  collector 
sPo  porosity 

o  characteristic  diameter  (A) 

p  total  mass  density  (kg  nrr3 ) 

Pane  density  of  anode  current  collector  (kg  m-3 ) 

Peat, ref  density  of  catalysts  which  actually  contribute  to  the 

reforming  reactions  (kgm-3) 

Pcat,sf  density  of  catalysts  which  actually  contribute  to  the 

water-gas  shift  reactions  (kg nrr3) 
r  twist  factor 

X  molar  fraction 

£2  collision  integration 

oji  mass  fraction  of  species  i 

Subscripts 

an  fuel  (anode  gas) 

ca  air  (cathode  gas). 


by  partially  reformed  methane  as  the  fuel.  From  the  compositional 
analysis  of  the  anode  gas,  their  experimental  results  indicated  that 
the  molar  ratio  of  the  hydrogen  and  carbon  elements,  i.e.,  H/C,  is 
not  a  constant  value  between  the  inlet  and  outlet  of  the  anode 
gas  channel,  and  that  nitrogen,  which  is  not  surely  present  in  the 
inlet,  is  significantly  detected  near  the  cell  edge.  The  following  con¬ 
clusions  are  suggested:  (1)  there  are  different  diffusion  velocities, 
or  different  species  velocities,  between  each  chemical  species  that 
composes  the  anode  gases,  and  (2)  reverse  diffusion  of  ambient  air 
from  the  periphery  of  the  cell  takes  place. 

In  this  paper,  to  theoretically  clarify  these  phenomena,  a  numer¬ 
ical  model  that  takes  the  convection  of  the  gas  and  diffusion  into 
account  was  constructed  for  a  seal-less  disk-type  SOFC  which  is 
operated  by  a  fuel  including  non-reformed  methane.  Therefore, 
we  compared  the  numerical  results  with  the  experimental  data 
such  as  the  concentration  distributions  of  the  gaseous  species 
and  voltage-current  characteristics  using  the  calculation  corre¬ 
sponding  to  the  experimental  conditions,  and  then  examined  these 
phenomena.  This  model  can  explain  the  influence  of  an  internal 
reforming  reaction  and  a  water-gas  shift  reaction  in  the  anode  gas 
flow  together  with  the  reverse  diffusion  from  the  periphery  of  the 
cell.  Furthermore,  this  model  suggests  the  distributions  of  an  elec¬ 
tromotive  force  and  a  current  because  it  deals  with  the  distributions 
of  the  concentration  profile  of  gaseous  species  along  the  flow  direc¬ 
tion,  thereby  the  voltage-current  characteristics  are  predicted. 

2.  Modeling 

2.1.  Simulation  model  and  assumption 

Fig.  1  shows  a  schematic  diagram  of  a  seal-less  disk-type  SOFC 
single  cell  as  a  simulation  model  used  to  explain  the  data  obtained 
from  the  experiment  of  Momma  et  al.  [20].  The  model  consists  of 
a  composite  disk  made  of  a  solid  electrolyte  with  electrodes  on 
both  sides,  and  current  collectors  and  metal  separators,  each  for 
the  positive  and  negative  electrodes. 
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I  metal  separator  I  cathode  current  collector  [ 

Fig.  1.  A  schematic  diagram  of  the  seal-less  disk-type  SOFC  single  cell. 


The  fuel,  which  consists  of  equilibrium  compositions  deter¬ 
mined  by  the  pre-reformer  temperature,  is  first  supplied  to  the 
anode  gas  channel  (current  collector)  through  orifices  located  at 
the  center  of  the  separator.  The  fuel  is  then  axial  symmetrically 
distributed  in  the  current  collectors  (Fig.  1 ).  The  air  is  also  supplied 
and  distributed  in  the  cathode  gas  channel  in  the  same  manner  as 
the  fuel.  The  remaining  fuel  after  the  fuel  cell  reaction  is  burned 
around  the  cell. 

The  following  describes  the  prerequisites  and  assumptions  of 
the  model: 


1.  The  gases  have  one-dimensional  axisymmetric  flows. 

2.  The  cell  temperature  is  constant  along  with  the  thickness  and 
radius  direction. 

3.  The  inner  pressure  of  the  cell  is  constant. 

4.  The  height  of  the  gas  channel  (the  current  collectors)  is  constant. 

5.  The  structure  of  the  cell  is  constant  throughout  the  cell. 

6.  The  composition  of  gases  consists  of  CH4,  H20,  H2,  CO,  C02  and 
N2  in  the  anode,  while  of  02  and  N2  in  the  cathode. 

7.  The  concentration  of  each  gas  species  is  constant  in  the  direction 
of  gas  channel  thickness. 

8.  As  an  inlet  boundary  condition,  at  the  center  of  the  disk-type 
cell,  the  composition  of  each  chemical  species  in  the  anode  gas 
is  in  equilibrium  as  determined  by  the  outlet  gas  temperature  of 
the  pre-reformer. 

9.  As  an  outlet  boundary  condition,  at  the  periphery  of  the  disk- 
type  cell,  the  air  exists  as  a  stoichiometric  amount  for  burning  all 
the  CH4,  H2,  and  CO  in  the  anode  gas.  In  other  words,  the  partial 
pressures  of  CH4,  H2  and  CO  are  zero,  and  N2  exists  0.79/0.21 
times  as  much  as  02  used  to  stoichiometrically  burn  the  fuel  at 
the  periphery  of  the  cell. 

For  items  1-5,  it  is  ideally  presumed  to  simplify  the  cell  model 
simulating  the  experiment.  As  for  item  1,  the  previous  experimen¬ 
tal  result  by  Momma  et  al.  [17]  significantly  demonstrated  that  the 
gaseous  flow  has  no  substantial  variations  along  the  circumferential 
direction  in  the  practical  seal-less  disk-type  cells  using  only  hydro¬ 
gen,  steam  and  nitrogen  as  the  anode  gases.  As  for  item  6,  oxygen  in 
the  anode  gas  can  be  ignored  on  the  basis  of  the  experiment.  Item 
7  is  an  assumption  for  the  one-dimensional  calculation,  meaning 
that  the  mixing  velocity  in  the  direction  of  the  thickness  is  much 
faster  than  the  velocity  of  the  radial  convection.  However,  there  is 
no  validity  near  the  inlet,  and  this  point  will  be  discussed  later.  The 
basis  of  item  8  is  a  result  of  the  experiment  [20].  Although  item  9 
is  a  bold  assumption  for  simplification,  it  is  confirmed  to  be  appro¬ 
priate,  because  the  simulation  result  and  experimental  result  well 
agree  as  will  be  shown  later. 

In  the  anode  current  collector,  the  fuel  containing  residual 
methane  and  carbon  monoxide  is  reformed  by  internal  reforming 


Table  1 

Description  of  the  anode  current  collectors  used  in  the  experiments. 


Type 

Material 

Thickness  [mm] 

Mass  [gnrr2] 

Density 
[kg  m~3] 

Type-A 

Ni-foam 

0.62 

440 

709 

Type-C 

Ni-SDC  impregnated 
Ni-foam/Ni-foam 

1.0 

1600 

1600 

reactions. 

Reforming  reaction:  CH4  +  H20  ->  3H2  +  CO  (1) 

Water-gas  shift  reaction:  CO  +  H20  ->  H2  +  C02  (2) 

The  reaction  rates  for  the  reforming  reaction  and  water-gas  shift 
reaction  are  assumed  in  Eqs.  (3)  and  (4),  respectively. 

>'ch4  =  1  -76  x  103exp  ^  A^ref  J  (pc h4  f'2  (3) 

where  rc h4  is  the  rate  of  the  reforming  reaction  [mols-1  kgCat— 1  ], 
A Earef  is  the  activation  energy  (=57,840J  mol-1 ),  R  is  the  gas  con¬ 
stant  (=8.314J  mol-1  K-1 ),  T  is  the  temperature  [K],  and  pCH4  is  the 
partial  pressure  of  methane  in  the  mixture  gas  [atm]  [21  ].  The  rea¬ 
son  why  Eq.  (3)  is  chosen  as  the  reaction  rates  for  the  reforming 
reaction  is  that  it  offers  a  good  fit  to  the  experimental  gas  concen¬ 
tration  distributions  as  compared  with  the  expression  reported  by 
Lee  et  al.  [22]. 

rco  =  26.1  x  103  exp  ^  (cco)11(ch2o)°'53(1  -  P)  (4) 


P  = 


Pco2Ph2 

Ks/PcoPh2o 


(5) 


where  rc o  is  the  rate  of  the  water-gas  shift  reaction 
[mol  s-1  kgcat-1  L  A EaSf  is  the  activation  energy  (=95,000J  mol-1 ), 
cCo  is  the  CO  concentration  in  the  mixture  gas  [molm-3],  cH2o 
is  the  H20  concentration  in  the  mixture  gas  [molm-3],  I<sf  is  the 
equilibrium  constant  of  Eq.  (2),  and  pz-  is  the  partial  pressure  of  the 
chemical  species  i  [23].  Therefore,  the  molar  reaction  quantity  per 
each  unit  time  and  unit  volume  in  Eqs.  (3)  and  (4)  [mols-1  m-3], 
cref  and  cSf  ,  are  obtained  from, 


Cref  ~  rCH4  Peat, ref  (6) 

CSf  =  rCO  Peat, sf  (7) 

where  pcat,ref  and  Pcat,sf  are  the  density  of  the  catalysts,  respec¬ 
tively  which  actually  contribute  to  the  reforming  and  water-gas 
shift  reactions  [kgm-3],  and  are  stated  as  follows: 


Peat, ref  —  Pane  &  ref  (&) 

Pcat,sf  —  Pane  OL sf  (9) 

where  pane  is  the  density  of  the  anode  current  collector  (see  Table  1 ) 
[kg  m-3  ],  and  aref  and  aSf  are  defined  here  as  the  utilizations  of  the 
catalyst  in  percentages  by  which  the  catalyst  in  the  anode  current 
collector  actually  contributes  to  the  above  reactions.  These  values 
are  experimentally  determined  as  shown  in  Section  3.1.1. 

On  the  other  hand,  the  concentrations  of  the  chemical  species  in 
the  fuel  and  air  change  based  on  the  electrochemical  cell  reactions 
in  Eq.  (10) 

H2  +  O2-  ->>  H20  +  2e-  at  anode  ,  . 

^02  +  2e-  O2-  at  cathode 

It  is  assumed  here  that  CO  is  consumed  only  in  the  shift  reaction 
and  not  by  the  direct  electrochemical  reaction  of  CO.  The  local  elec¬ 
tromotive  force  Eemf  is  obtained  by  the  Nernst  equation  from  the  02 
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Fig.  2.  Equivalent  circuit  of  the  disk-type  cell. 


partial  pressure  ratio  of  the  anode  and  cathode  assuming  that  H2, 
02  and  H20  in  the  anode  gas  are  in  thermal  equilibrium. 


Eemf  ~ 


a  g; 


h20 


RT  , 


P^Pof 


2  F 


+  ^ln 
2  F  pH?o 


(11) 


where  AG^q  is  the  Gibbs  free  energy  (H2  +  02  =  H20),  and  Fis  Fara¬ 
day  constant. 

The  equivalent  circuit  of  the  disk-type  cell  assumed  as  the  model 
is  shown  in  Fig.  2  [17],  which  stands  for  a  local  circuit  segmented 
from  the  whole  SOFC.  Each  symbol  in  Fig.  2  denotes  as  follows:  Rif 
ohmic  resistance  of  the  electrolyte  by  ionic  conduction  [to2];  Re , 
ohmic  resistance  of  electrolyte  by  electronic  conduction  [£2m2];  Rr , 
reactive  resistance  of  the  electrodes  [£2m2];  Rs,  series  resistance 
of  the  electrodes  [to2];  f,  ionic  current  density  which  flows  in 
the  electrolyte  [Am-2];  Je,  electronic  current  density  which  flows 
in  the  electrolyte  [Am-2];  J0,  current  density  which  flows  in  an 
external  circuit  [Am-2];  and  Vceu,  cell  operating  voltage  [V].  These 
cell  resistances  are  assumed  to  be  a  function  of  the  temperature, 
and  are  given  by 


Rcell  ~  Pcell0 


~Ear 

(-- 

1  'll 

~Rf 

It 

1073  7. 

(12) 


where  Rcen  means  each  resistances  described  above  in  [£2m2],  Rcem 
is  the  resistance  at  1073  K  [£2m2],  and  Ear  is  the  activation  energy 
[J  mol-1  ]  (see  Table  2)  [8].  These  values  were  experimentally  deter¬ 
mined  by  an  AC  impedance  measurement  of  the  cell.  The  whole  cell 
current  is  written  as 

rrcell 

Cell  =  /  Jo2jrr  dr  (13) 

Jo 

where  J0  means  the  local  current  density  and  is  obtained  from 


J  _  EemfP-e  -  Vcell(Pi  +  Pr  +  Re)  .  . 

J°  ~  Re(Ri  +Rr+Rs)  +  Rs(Ri  +Rr)'  1  j 

The  molar  reaction  rate  in  the  anode  and  cathode  by  the  cell 
reaction  for  each  unit  time  and  unit  volume  [mols-1  m-3],  Cj  and 
cj,  is  obtained  from,  respectively 


•  _Jo+Je  1  _  Ji  1 

J  2F  dan  2  Fdan 

Jo+Je  1  _  Jj  1 

j  2  F  dca  2Fdca 


(15) 

(16) 


where  dan  is  the  thickness  of  the  anode  current  collector  (see 
Table  2)  and  dca  is  the  thickness  of  the  cathode  current  collector. 


Table  2 

Simulation  parameters. 
Cell  size 


Radius  of  the  cell 

rcell 

6  cm 

Thickness  of  the  anode  current  collector 

dan 

See  Table  1 

Porosity  of  the  anode  current  collector 

&an 

0.95 

Thickness  of  the  cathode  current 

dca 

0.07  cm 

collector 

Porosity  of  the  cathode  current  collector 

&ca 

0.95 

Physical  properties  of  cell 

Density  of  the  anode  current  collector 

Pan 

See  Table  1 

Twist  factor  of  the  anode  current 

tan 

1 

collector 

Twist  factor  of  the  cathode  current 

~Cca 

1 

collector 

Ohmic  resistance  of  electrolyte  by  ionic 

Rw 

0.17  Q  cm2 

conduction  at  800  °C 

Activation  energy  of  R; 

Eai 

54,200 J  mol”1 

Ohmic  resistance  of  electrolyte  by 

ReO 

20  Q  cm2 

electron  conduction  at  800  °C 

Activation  energy  of  Re 

Eae 

54, 200 J  mol”1 

Reactive  resistance  of  the  electrodes  at 

RrO 

0.04  £2  cm2 

800  °C 

Activation  energy  of  Rr 

Ear 

50, 200 J  mol-1 

Series  resistance  of  the  electrodes  at 

Rso 

0.002  £2  cm2 

800°C 

Activation  energy  of  Rs 

Eas 

54, 200 J  mol-1 

Operating  conditions 

Flow  rate  of  the  fuel  (methane) 

85  Nccmin-1 

Flow  rate  of  the  air 

1700  Ncc  min-1 

Oxygen  fraction  in  the  air 

0.21 

Inner  pressure  of  the  cell 

P 

1  atm 

Room  temperature 

Troom 

25  °C 

Outlet  gas  temperature  of  pre-reformer 

Tref 

650,  600,  550 °C 

Temperature  of  the  cell 

T'cell 

750°C 

Fuel  utilization 

Vf 

0,  50,  70% 

Steam-carbon  ratio 

S/C 

3.0 

2.2.  Transport  equation  of  gas  flow 

The  concentration  distribution  of  each  gaseous  species,  cell  volt¬ 
age  and  cell  current  density  are  calculated  in  such  a  way  that  the 
equations  of  mass  balance  and  the  equation  of  continuity  are  con¬ 
verted  into  discrete  equations  by  the  control-volume  method  using 
the  implicit  method  and  upwind  scheme. 

2.2.1.  Equation  of  mass  balance  for  each  species 
Anode  side: 

CH4  :  —  Pan ^CH4 tlan  ~  EDCH4,mPan  ^ 4  ^  =  (~Cref)  Mc h4 

(17) 


i_r  n  .  9  /  n  da)H2o  \ 

H2O  .  rfjp  I  r PanOJ^Otlan  ~  rDH20,mPan — gp —  j 

=  (~Cref  ~  Csf  +  9)  ^H20  (18) 


H2  : 


d  f  dcou2 

jrjp  I  r panO)\\2uan  ~  r^H2 , m Pan  ^ 


=  (3 Cref  +  Csf  -  Cj)  Mh2 


(19) 


CO  : 


_d_ 

rdr 


panCOcoUan  ~  tDc o  ,mPan 


dajco 

~dr~ 


(cref  -  csf)  Mco 

(20) 


12 


T.  Shimada  et  al.  /  Journal  of  Power  Sources  187  (2009)  8-18 


9  (  dcoco-y  \  / .  \ 

C02  :  —  I  r pan(t>c02uan  -  tDqq2  m Pan  I  =  [Csf )  ^C02 

(21) 

N2an  •  PanCON2anUan  ~  r^N2an,mPan — =  0  (22) 

Cathode  side: 

9  (  9&>o9  \  /  \ 

°2  *  ^  (  rpca(Oo2uca  -  rDo2N2Pca-^  J  =  (— 0.5cj)  Mq2  (23) 

9  /  9&)NoCa\ 

^2ca  •  (  rPca&>N2cauca  —  r^N202  Pea — ^ —  1=0  (24) 

where  p  is  the  total  mass  density  [kgnrr3],  is  the  mass  fraction 
of  species  i  [-],  u  is  the  mass  average  velocity  [ms-1],  Dim  is  the 
diffusion  coefficient  of  species  i  in  the  multi-component  mixture 
gas  [m2s-1],  Mt  is  the  molecular  weight  of  species  i  [kg  mol-1], 
suffix  an  means  the  fuel  (anode  gas)  and  suffix  ca  the  air  (cathode 
gas). 


where  DAB  is  the  mutual  diffusion  coefficient  of  species  A  and 
B  [m2  s-1  ],  Ma  and  MB  are  the  molecular  weights  of  species 
A  and  B  [gmol-1],  respectively,  P  is  the  pressure  [atm],  oA 
and  oB  are  the  characteristic  diameters  of  species  A  and  B 
[A],  respectively,  £2D  is  the  collision  integration,  T*  =  kT/sAB, 
eAB  =  (eA  £b)1/2,  >1  =  1.06036,  B  =  0.15610,  C=0.19300,  D  =  0.47635, 
£=1.03587,  £=1.52996,  G=  1.76474,  H  =  3.89411,  e/k  is  the  inter- 
molecular  force  constant  [K]  and  k  is  Boltzmann  constant 
(=1.3807  xlO-23JI<-1)  [25]. 

An  effective  diffusion  coefficient,  De,  in  porous  current  collectors 
is  calculated  as  follows: 

De  =  ~~Dab  (34) 

where  r  is  the  twist  factor  [-]  and  £po  is  the  porosity  [-]  [26].  In 
this  paper,  the  equations  of  mass  balance  shown  in  Eqs.  (17)-(24) 
are  calculated  by  taking  the  effective  diffusion  coefficient  into  con¬ 
sideration. 

2.3.  Analytical  conditions 


2.2.2.  Equation  of  continuity 

The  equation  of  mass  continuity  is  as  follows.  Derivation  of  the 
expression  is  described  in  Appendix  A: 


IT  +  Wr{rPanUan)  - 

(1cj)m o2  at  anode 

(25) 

d~W  +  Wr[rpcaUca)=\ 

(-1c;)  Mq2  at  cathode 

(26) 

The  mass  density  of  the  fuel  and  air  is  given  by 

Pan  —  J2CanXM 
i 

(27) 

Pea  =  'y^CcaXjMj 

(28) 

j 


where  x  is  the  molar  fraction  [-],  and  the  molar  density  [mol  m-3  ], 
c,  is  a  value  dependent  only  on  the  temperature  T  and  pressure  P, 
which  are  assumed  to  be  constant  in  the  cell. 


(29) 


2.2.3.  Diffusion  coefficient 

The  diffusion  coefficients  of  a  chemical  species  in  the  anode  and 
cathode  used  in  Eqs.  (17)-(24)  are  calculated  from  the  prediction 
value  in  a  multi-component  gas  compound  [24]. 


1  -X/  =  *j 
Hi  m  'Hjj 

j= 2 


(30) 


where  D lm  is  the  diffusion  coefficient  when  element  1  diffuses  in 
a  mixture  gas  m,  and  Dy  is  the  diffusion  coefficient  of  the  pair  1  -j 
of  species  in  a  binary  mixture. 

The  diffusion  coefficient  of  a  binary  mixture  gas  composed  of 
A  and  B  at  low  (normal)  pressure  is  calculated  by  the  following 
expressions: 


Dab  =  1.858  x  10 
oA  +  aB 

&AB  =  - n - 


— 3T3/2  [(m^  +  mb)/mamb]1/2 

P<nD 


&D  —  + 


+ 


+  ' 


T*b  exp  (DT*)  exp(FT*)  exp(HT*) 


(31) 

(32) 

(33) 


In  the  experiment  by  Momma  et  al.,  three  types  of  current  col¬ 
lectors  with  different  structures  and  materials,  which  also  serve  as 
the  channel  for  the  fuel  gas,  were  tested  for  comparison  of  the  inter¬ 
nal  reforming  characteristics  [20].  However,  only  two  types  of  the 
current  collectors  having  significantly  different  reforming  charac¬ 
teristic,  type-A  and  type-C  (see  Table  1 ),  are  discussed  in  this  paper 
because  the  reforming  characteristic  of  the  current  collector  type-A 
was  similar  to  that  of  type-B.  Measurements  were  made  by  vary¬ 
ing  the  fuel  utilization  (Uf=  0,  50,  70%)  at  various  pre-reformer  exit 
temperatures  (Tref=  550,  600,  650  °C)  with  a  steam  to  carbon  ratio 
of  S/C  =  3.0  at  the  cell  temperature  of  750  °C. 

The  physical  property  values  and  the  operating  conditions  used 
in  the  calculations  are  indicated  in  Table  2.  The  simulation  param¬ 
eters,  such  as  cell  sizes  and  physical  properties  and  operating 
conditions,  are  adjusted  to  the  experimental  conditions  as  shown 
in  reference  [20]. 

3.  Results  and  discussion 

33.  Anode  current  collector  type-C 

333.  Decision  on  parameters  of  the  reaction  rate  by  comparison 
with  experimental  results 

As  for  the  reaction  rates  of  the  reforming  reaction  and  water- 
gas  shift  reaction,  the  percentages  that  the  catalyst  (Ni)  actually 
works  for  the  reactions  are  treated  as  the  unknown  parameters, 
aref  and  asf.  The  optimum  values  to  fit  the  calculated  gas  concen¬ 
tration  profile  to  the  experimental  ones  are  found  by  varying  aref 
and  aSf  at  only  one  condition  of  Tref=  550  °C  and  Uf=  70%.  The  reason 
why  Tref=  550  °C  is  chosen  is  that  more  non-reformed  methane  is 
contained  in  the  fuel  cell  due  to  the  reforming  temperature  lower 
than  650  and  600 °C.  Also,  the  reason  why  Uf  =70%  is  chosen  is 
that  the  greatest  variation  in  each  chemical  species  concentration 
is  observed  between  the  center  and  the  periphery  of  the  cell  due 
to  the  most  progressive  electrochemical  reactions.  Therefore,  it  is 
an  advantage  to  confirm  the  effect  of  the  selected  parameters  on 
the  gas  concentration  distributions  when  the  numerical  results  are 
compared  to  the  experimental  data.  The  parameters  for  the  per¬ 
centage  and  the  densities  corresponding  to  it  are  listed  in  Table  3. 
The  value,  aref,  is  assumed  to  be  equal  to  aSf  for  simplification  of 
the  calculation. 

The  gas  concentration  distributions  of  CH4,  H20,  H2,  CO  and  C02 
for  each  condition  in  Table  3  are  shown  in  Figs.  3-7.  The  experi¬ 
mental  values  at  the  radii  of  0,  1.4,  2.8,  4.2  and  5.6  cm  are  jointly 
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Table  3 

Conditions  for  a  decision  on  the  reaction  rate  parameters. 


Radius  distance  /  cm 

Fig.  3.  CH4  distributions  for  each  condition  and  experimental  data  (Tref=  550  °C, 
Uf=  70%,  type-C). 

plotted  on  the  figures.  It  can  be  seen  from  these  figures  that  the 
experimental  values  are  the  closest  to  the  calculation  results  for 
condition  3,  i.e.,  aref= asf=  0.03,  or  pre/=Ps/=48 kgm-3.  Therefore, 
the  density  of  the  catalysts,  which  actually  contribute  to  the  reform¬ 
ing  and  water-gas  shift  reactions  is  set  at  48  kg  m-3  for  the  anode 
current  collector  type-C.  Using  this  value,  the  gas  concentration  dis¬ 
tributions  are  likewise  calculated  for  the  other  conditions,  and  the 
validity  of  the  parameter  is  evaluated  in  the  following  section. 


Radius  distance  /  cm 

Fig.  4.  H20  distributions  for  each  condition  and  experimental  data  (Tref=  550  °C, 
Uf  =70%,  type-C). 


Radius  distance  /  cm 

Fig.  5.  H2  distributions  for  each  condition  and  experimental  data  (Tref=  550  °C, 
Uf =70%,  type-C). 

As  a  further  consideration,  from  Fig.  3,  the  experimental  CH4 
concentration  hardly  changes  between  r  =  0  cm  and  r  =  1.4  cm  unlike 
the  calculated  result,  and  the  experimental  H20  concentration  at 
1.4  cm  is  higher  than  the  calculated  one,  while  the  experimen¬ 
tal  H2  and  CO  concentrations  are  lower  than  the  simultaneously 
calculated  values  (see  Figs.  4-6).  It  is  interpreted  that  there  is  no 
reforming  reaction  in  this  region.  That  is,  the  region  where  the  cat¬ 
alyst  effectively  performs  becomes  small  because  the  diffusion  of 
the  gas  in  the  thickness  direction  is  lower  than  the  radial  convective 
flow  due  to  the  higher  velocity  near  the  center.  On  the  other  hand, 
in  this  model,  the  diffusion  in  the  thickness  direction  is  assumed  to 
be  sufficiently  greater  than  the  radial  convection.  In  other  words, 


Radius  distance  /  cm 

Fig.  6.  CO  distributions  for  each  condition  and  experimental  data  (Tref=  550  °C, 
1^=70%,  type-C). 
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Fig.  7.  CO2  distributions  for  each  condition  and  experimental  data  (Tref=  550  °C, 
Uf=  70%,  type-C). 

the  assumption  that  the  distribution  of  the  gas  concentration  in 
the  thickness  direction  is  not  considered  is  not  suitable  in  the  near 
inlet. 

However,  the  entire  tendency  of  the  numerical  results  is  in  good 
agreement  with  the  experimental  values  although  the  distribution 
is  not  in  close  agreement  from  r  =  0  to  r=  1.4  cm. 

3A.2.  Gas  concentration  distributions  in  the  cell  at  various 
pre-reformer  temperatures  and  fuel  utilization 

Using  the  reaction  rate  parameter  determined  in  the  above 
investigation,  the  distributions  of  gas  concentrations  are  calculated 
and  then  compared  to  the  experimental  results  for  each  condition  of 


Fig.  8.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and  cur¬ 
rent  density  ( Tref=  650  °C,  Uf=  0%,  type-C). 


Radius  distance  /  cm 

Fig.  9.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and  cur¬ 
rent  density  ( Tref  =  650  °C,  Uf=  50%,  type-C). 

the  pre-reformer  temperatures  of  650, 600  and  550  °C  and  fuel  uti¬ 
lizations  of  0,  50  and  70%.  Thereafter,  the  validity  of  the  simulation 
model  is  examined.  The  distributions  of  each  gas  concentration, 
electromotive  force,  cell  voltage  and  current  density  are  displayed 
in  Figs.  8-15.  The  result  when  Tref=  550  °C  and  Uf=  50%  is  omitted 
because  it  was  not  measured.  The  experimental  values  at  the  radii 
of  0,  1.4,  2.8,  4.2  and  5.6  cm  are  jointly  plotted  on  the  figures  indi¬ 
cated  by  the  symbols  of  (O)  CH4,  (□)  H20,  (A)  H2,  (v)  CO,  (0)  C02 
and  (■)  N2  in  the  anode. 

The  calculated  concentration  distributions  of  the  anode  gas 
species  well  agree  with  the  experimental  in  every  case  for  the 
pre-reformer  temperatures  and  the  fuel  utilizations. 


Fig.  10.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  ( Tref=  650  °C,  Uf=  70%,  type-C). 
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Fig.  11.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  (Tref  =  600  °C,  Uf=  0%,  type-C). 


Radius  distance  /  cm 

Fig.  13.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  ( Tref=  600  °C,  Uf=  70%,  type-C). 


In  the  experiment,  the  distinct  effect  of  internal  reforming  is 
observed  for  the  anode  current  collector  type-C.  As  shown  in  every 
figure,  the  experimental  well  agrees  with  the  calculation  results  of 
the  CH4  concentration  that  decreases  along  the  flow  direction  due 
to  the  reforming  reaction.  This  good  agreement  suggests  that  the 
reforming  reaction  of  the  cell  used  in  the  experiment  can  be  well 
described  by  the  reaction  rate  equation  (Eq.  (3)). 

For  the  CO  and  C02  concentration  distributions,  the  considered 
model  can  describe  the  phenomenon  that  the  CO  concentration 
increases  (C02  decreases),  i.e.,  the  reverse-shift  reaction,  from  the 
center  of  the  cell  to  around  1.4  cm  in  radius  (2.8  cm  in  Uf=  0%  due  to 
a  lower  H20  concentration),  while  the  CO  concentration  decreases 


Radius  distance  /  cm 

Fig.  12.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  (Tref  =  600  °C,  Uf=  50%,  type-C). 


(C02  increases),  i.e.,  the  shift  reaction,  from  its  outside  because  the 
H20  concentration  increases  and  the  H2  concentration  decreases 
due  to  the  electrochemical  reaction.  Therefore,  the  reaction  rate 
equation  of  the  shift  reaction  can  model  the  experimental  result. 

It  can  be  confirmed  from  the  calculation  results  that  the  N2  con¬ 
centration  distributions  observed  in  the  anode  gas  channel  in  the 
experiment  are  due  to  the  reverse  diffusion  from  the  periphery.  It  is 
noted  from  Figs.  8  to  15  that  the  reverse  diffusion  of  N2  takes  place 
more  remarkably  as  the  fuel  utilization  Uf  is  lower;  i.e.,  hydrogen 
exists  more  in  the  gas  channel  due  to  its  lower  consumption.  This  is 
because  the  diffusion  coefficient  in  the  anode  gas  mixture  of  nitro¬ 
gen  becomes  high  when  the  fraction  of  the  hydrogen  component 
in  the  gas  increases.  This  is  also  further  discussed  in  Section  3.1.3. 


Fig.  14.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  ( Tref=  550  °C,  Uf= 0%,  type-C). 
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Fig.  15.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  (Tre/=  550  °C,  Uf=  70%,  type-C). 

When  Uf=  0%,  the  concentration  of  H2  sharply  decreases  in  the 
range  r>4.2cm.  This  is  caused  by  the  decrease  in  the  H2  concen¬ 
tration  at  the  outlet  cell  because  H2  is  consumed  by  the  oxidization 
with  02  in  the  periphery  of  the  cell,  leading  to  an  increased  H2 
diffusion  velocity  by  the  concentration  difference.  The  boundary 
condition  at  the  cell  outlet  is  assumed  for  H2  to  be  completely  con¬ 
sumed.  The  calculated  H2  and  N2  concentrations  at  r  =  5.6  cm  are  in 
good  agreement  with  the  experimental  data  for  each  fuel  utiliza¬ 
tion.  Therefore,  the  assumption  at  the  outlet  boundary  condition 
can  describe  the  experiment  results  very  well.  As  a  consequence, 
as  shown  in  Figs.  8-15  it  is  predicted  that  the  electrolysis  current  is 
induced  in  the  vicinity  of  the  cell  edge  in  which  Eemf  becomes  less 
than  Vceu. 

Therefore,  it  is  summarized  that  the  reverse  diffusion  of  air  from 
the  cell  edge,  the  internal  reforming  reaction  and  the  water-gas 


shift  reaction  can  be  well  described,  resulting  in  a  good  agreement 
between  the  calculated  and  experimental  results  for  the  concen¬ 
tration  distributions  of  the  anode  gases. 

3  A3.  Calculation  ofH/C  mass  balance,  and  comparison  to 
experiment  results 

The  carbon  balance,  Cbalance,  and  hydrogen  mass  balance, 
^balance*  are  still  disputed.  Cbalance,  and  Hbalance,  are  defined  by  Eq. 
(35)  and  Eq.  (36),  respectively. 

^balance  =  XCH4  +  XCO  +  XC02  (35) 

^balance  =  4XCH4  +  2Xh2  +  2Xh20  (36) 

where  Xi  denotes  the  molar  fraction  of  gaseous  species  i. 

In  Fig.  16,  the  calculated  Hbaiance /Cbalance  is  plotted  together  with 
the  experimental  values.  The  calculation  results  show  a  tendency 
similar  to  the  experimental  values  versus  the  various  parame¬ 
ters,  such  as  r,  Tref  and  Uf.  For  the  experimental  values,  while  the 
observed  Hbalance /Cbalance  remains  approximately  equal  to  its  input 
value  of  10  at  4.7  cm,  it  suddenly  decreases  at  r  =  56cm.  Addi¬ 
tionally,  it  decreases  more  with  lower  fuel  utilization.  This  can 
be  explained  as  follows.  The  low  fuel  utilization  results  in  the 
significant  existence  of  hydrogen  which  causes  a  high  diffusion 
coefficient  of  hydrogen  and  high  concentration  gradient,  dcou2/dr , 
in  the  vicinity  of  the  cell  edge,  resulting  in  the  faster  molecular  aver¬ 
age  velocity  of  hydrogen  gas  than  the  mass  average  velocity.  On  the 
other  hand,  the  molecular  average  velocity  of  the  other  chemical 
species  including  carbon  is  much  slower  than  that  of  hydrogen  gas. 
The  concentration  of  the  chemical  species  increases  when  its  veloc¬ 
ity  decreases  because  of  the  continuity  law,  and  the  concentration 
decreases  when  its  velocity  increases.  The  ratio  of  FIbalance/Cbaiance 
then  decreases  more  at  r  =  56  cm  with  the  lower  fuel  utilization. 

Moreover,  for  the  same  fuel  utilization,  the  Hbalance / Cbalance  is 
greater  as  the  exit  temperature  of  the  pre-reformer  is  lowered 
because  the  non-reformed  methane  increases,  thus  resulting  in  a 
hydrogen  decrease. 

As  shown  above,  the  experimental  result  that  the  ratio  of  H/C 
does  not  remain  a  constant  value  in  the  anode  gas  for  the  seal¬ 
less  disk-type  SOFC  can  be  explained  by  a  numerical  model  that 
considers  both  the  convection  and  diffusion  of  the  gas. 
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Fig.  16.  H/C  mass  balance. 


Fig.  17.  Voltage-current  characteristics. 
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Radius  distance  /  cm 

Fig.  18.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  (Tref  =  650  °C,  Uf= 0%,  type-A). 

3.1.4.  Voltage-current  characteristics 

Fig.  17  shows  a  comparison  of  the  voltage-current  (V-/)  char¬ 
acteristic  between  the  calculation  and  the  experiment  results.  In 
the  0-70%  range  of  fuel  utilization,  the  calculated  V-I  curves  are 
in  reasonably  good  agreement  with  the  experimental  plots  at  each 
temperature.  The  validity  of  the  simulation  model  can  also  be  veri¬ 
fied  for  the  cell  generating  performance  characteristic  as  well  as  the 
gaseous  chemical  reactions.  At  the  fuel  utilization  of  80%,  however, 
there  is  the  difference  from  the  experimental  plots. 

Also,  it  was  not  possible  to  calculate  up  to  a  90%  fuel  utiliza¬ 
tion  by  this  model  at  Tref=  650  and  600  °C  as  mentioned  above.  One 
likely  reason  for  this  difference  is  that  it  relates  to  the  electrolysis 
current  in  the  vicinity  of  the  outer  cell,  i.e.,  it  might  be  not  appre¬ 


Radius  distance  /  cm 

Fig.  19.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  (Tref=  650  °C,  Uf=  70%,  type-A). 


ciable  in  this  model  whether  corresponding  electrolysis  current  is 
obtained  between  the  calculated  and  experimental  results  or  not, 
because  of  writing  the  each  cell  resistance  by  the  expression  of  the 
simple  temperature  dependence.  Another  is  that  the  electrochem¬ 
ical  direct  reaction  of  CO  might  have  taken  place  during  the  actual 
reaction,  although  the  reaction  is  assumed  not  to  in  this  model. 
This  is  also  further  discussed  in  Section  3.2.  The  development  of  a 
precise  model  about  these  points  is  future  work. 

3.2.  Anode  current  collector  type-A 

Similar  to  the  anode  current  collector  type-C,  aref  and  asf 
are  determined  by  fitting  the  distributions  of  the  gas  concentra¬ 
tion  at  Tref=  600  °C  and  Uf  =70%.  As  the  sequence  for  comparison 
with  the  experimental  data,  the  following  catalytic  parameters 
were  obtained;  aref=  0.03  and  aSf=  0.005,  i.e.,  pref=  21.2  kg  m-3  and 
pSf=  3.55  kg  m-3.  This  corresponds  to  the  effective  amount  of  cata¬ 
lyst  of  12.7  g  m-2  for  the  reforming  reaction  while  that  of  type-C  is 
48  g  m-2.  Although  an  internal  reforming  hardly  occurs  in  the  type- 
A  during  the  experiment,  it  was  found  that  the  amount  of  catalyst 
per  unit  area,  which  contributes  to  the  internal  reforming  reaction, 
is  almost  1  / 4  of  that  for  the  type-C. 

By  using  these  values,  the  gas  concentration  distributions  are 
calculated  and  then  compared  with  the  experimental  results  for 
each  condition  when  the  pre-reformer  temperature  is  650  and 
600  °C  and  the  fuel  utilization  is  0,  50  and  70%.  The  distributions  of 
each  gas  concentration,  electromotive  force,  cell  voltage  and  cur¬ 
rent  density  when  Uf=  0,  70%  as  an  illustration  are  displayed  in 
Figs.  18-21.  From  the  figures,  for  the  type-A,  the  numerical  results 
are  also  in  good  agreement  with  the  experimental  for  every  condi¬ 
tions.  The  results  when  Uf=  50%  are  omitted  since  they  also  agree 
well  with  the  experimental  data  as  well  as  those  when  Uf=  0,  70%. 

When  Uf=  70%,  it  is  noted  that  the  calculated  CO  concentration 
is  higher,  while  the  calculated  C02  concentration  is  lower  than  the 
experimental  data  at  r  =  5.6cm.  It  is  suggested  that  the  electro¬ 
chemical  direct  reaction  of  CO  on  the  anode  might  have  taken  place 
during  the  actual  reaction  because  the  H2  concentration  is  signifi¬ 
cantly  decreasing  at  r  ^  4.2  cm  due  to  the  poor  performance  of  the 
internal  reforming. 


Radius  distance  /  cm 

Fig.  20.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  (Tref=  600  °C,  Uf= 0%,  type-A). 
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Fig.  21.  Distributions  of  gas  concentration,  electromotive  force,  cell  voltage,  and 
current  density  ( Tref=  600  °C,  Uf=  70%,  type-A). 


4.  Conclusions 

The  characteristics  of  the  electrical  power  generation  and  the 
internal  reforming  of  a  seal-less  disk-type  SOFC  cell  have  been 
numerically  analysed  by  a  simulation  to  explain  the  previously 
reported  experimental  results. 

From  comparisons  between  the  calculations  and  experimental 
results  for  the  anode  gas  concentration  profile  and  the  voltage- 
current  characteristics,  the  validity  of  the  simulation  model  has 
been  verified,  and  the  model  was  found  to  explain  the  influence 
of  the  internal  reforming  reaction  with  the  water-gas  shift  reac¬ 
tion,  the  reverse  diffusion  of  air  from  the  cell  circumference  and 
variations  in  the  molar  ratio  of  hydrogen  and  carbon  in  the  cell  out¬ 
let  region  for  a  disk-type  SOFC  cell  which  is  operated  by  partially 
reformed  methane. 


Appendix  A 

Eq.  (25)  is  derived  as  follows.  The  equation  of  mass  continuity 
in  the  anode  is  given  by 

d  (  \ 

foPan  CH4  +  &>H20  +  &>H2  +  wCO  +  &>C02  +  &>N2an  j 


9 

+  ^  [rPan  Kh4  +  ^H20  +  ^H2  +  &>CO  +  &>C02  +  &>N2an) 

doJcu4 


D{ 


-rpan 


CH 4m  ~ 


dr 


+  Dh 


do)R20  n  d(L>H2 
'H20 m  dr  T  L)H2m 

,  n  dcoc o  n  dcoCo2 

+UcO m^z - f  T'CO 2m 


dr 


dr 


■  +  ^2ai 


9r 

9^N2ai 


dr 


The  summation  of  the  mass  fraction  of  each  gas  and  that  of  the  mass 
movement  due  to  mutual  diffusion  are  written  as  follows: 

(&>ch4  +  &>h2o  +  &>h2  +  <*>co  +  &>co2  +  <^N2an )  =  1  (38) 


9^ch4  9^h20  dcou2 

UCH  4m  gr  +  ^H20  m  ;]r  +  T>H2m 


+  DCOm  —r^~  +  DCo2m 


dr 

dcoco2 

dr 


'2  anm 


dr 

9®N2ai 


dr 


=  0 


(39) 


Therefore,  the  equation  of  mass  continuity  shown  as  Eq.  (25)  is 
obtained,  substituting  Eqs.  (38)  and  (39)  into  Eq.  (37). 
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